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Abstract 

As the most important discovery channel for a light Higgs boson at the LHC, the di-photon 
signal gg > h > yy is sensitive to underlying physics. In this work we investigate such a signal 
in a comparative way by considering three different supersymmetric models, namely the minimal 
supersymmetric standard model (MSSM), the next-to-minimal supersymmetric standard model 
(NMSSM) and the nearly minimal supersymmetric standard model (nMSSM). Under the current 
collider and cosmological constraints we scan over the parameter space and obtain the following 
observation in the allowed parameter space: (i) In the nMSSM the signal rate is always suppressed; 
(ii) In the MSSM the signal rate is suppressed in most cases, but in a tiny corner of the parameter 
space it can be enhanced (maximally by a factor of 2); (iii) In the NMSSM the signal rate can be 
enhanced or suppressed depending on the parameter space, and the enhancement factor can be as 


large as 7. 


PACS numbers: 


I. INTRODUCTION 


So far the most important question to be answered in particle physics is the mechanism 
of the electroweak symmetry breaking and thus hunting for the Higgs boson responsible 
for it is the main task of current collider experiments. In the framework of the Standard 
Model (SM), the mass of the Higgs boson is preferred to be 116.4*159 GeV by precision 
electroweak data |l]. To search for such a relatively light Higgs boson, great efforts have 
been made in LEP and Tevatron experiments, which reported null results and excluded a 
Higgs boson with m; < 114.4 GeV and 158 GeV € mp < 175 GeV at 95% C.L.. 
The Large Hadron Collider (LHC) is more powerful in discovering the SM Higgs boson, and 
depending on its mass, different search strategies will be applied. For a light Higgs boson 
below 140 GeV, although its largest signal at the LHC is bb from the gluon-fusion process 
gg — h— bb li such a signal is undetectable due to the overwhelming QCD background; 
instead, the rare decay mode h > yy with Br(h > yy) œ 0.296 for mp, =120 GeV offers 
a very clean signature to make the di-photon signal gg — h — yy a promising discovery 
channel. It is now expected that, with 2fb-! integrated luminosity at the LHC running at 
Vs = 7 TeV, the di-photon signal is able to exclude the light Higgs boson in the SM [5]. 

In low energy supersymmetric models (SUSY), the SM-like Higgs boson (the CP-even 
Higgs boson with largest coupling to gauge bosons) is usually predicted with mass below 
about 140 GeV. For such a Higgs boson, although there may exist other discovery channels 
at the LHC, the di-photon channel gg — h — ^y is still one of the most important discovery 
modes. So, studying this signal will allow for a probe of low energy SUSY and, as emphasized 
in |6], even a discrimination of different models. Although in the literature some studies of 
the signal have been presented in SUSY , these analyses were performed separately in 
different models and a comparative study is necessary in order to discriminate the models. 
In this work we perform such a comparative study by considering three different SUSY 
models, namely the minimal supersymmetric standard model (MSSM), the next-to-minimal 
supersymmetric standard model (NMSSM) [10,|11] and the nearly minimal supersymmetric 
standard model (nMSSM) |12, d We will scan over the parameter space under current 
constraints from collider experiments and the neutralino dark matter relic density, and then 
in the allowed parameter space we calculate the di-photon signal rate and compare the 


results for different models. 


This work is organized as follows. We first briefly describe the three supersymmetric 
models in Sec. II. Then we present our numerical results and discussions in Sec. III. Finally, 


we draw our conclusions in Sec. IV. 


Il. SUPERSYMMETRIC MODELS 


As the most economical realization of SUSY in particle physics, the MSSM has been 
intensively studied. However, since this model suffers from the jj—problem and the little 
hierarchy problem, some of its extensions like the NMSSM and nMSSM were recently paid 


attention to[10]. The differences of these models come from their superpotentials: 


Waussm = Wr + pH, - Ha, (1) 
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Wyoussm = Wr + àH,- H4S + 355^. (2) 

Wamssm = Wr AH, - H4$ + £p M26, (3) 


where W is the MSSM superpotential without the u term, Fa and $ are the Higgs doublet 
and singlet superfields respectively, and the dimensionless coefficients A, & and € and the 
dimensional coefficients u and Mp are usually treated as independent parameters. In the 
NMSSM and nMSSM, when the scalar component (S) of the singlet Higgs superfield $ 
develops a vacuum expectation value (VEV), an desired effective -term (peg; = A(S)) is 
generated at the weak scale. Note that the nMSSM differs from the NMSSM in the last 
term of the superpotential, where the cubic singlet term S? in the NMSSM is replaced by 
the tadpole term €-M?. Considering that the tadpole term does not induce any interaction, 
one can infer that, except for the minimization conditions of the Higgs potential and the 
mass matrices of the Higgs bosons, the nMSSM is actually identical to the NMSSM with 
vanishing &. 


Corresponding to Eq. (L413), the soft-breaking terms in Higgs sector are given by 


Vee = mH + mH? + (Bul: Hat hc.), (4) 
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where Mu, Ma, Ms, B, A, and A, are all soft parameters. Like the usual treatment of the 


multiple Higgs theory, one can write the scalar fields H,, Ha and S as 


H+ tat eared 1 
A= i , Ha = k , S= — (s +0 +it), 7 
testes a H; V2 ( £) (7) 


and diagonalize the mass matrices of the Higgs bosons to get their mass eigenstates: 
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In above expressions, ^, ho, h3 and a, A denote physical CP-even and CP-odd neutral Higgs 
bosons respectively, G? and G^ are Goldstone bosons eaten by Z and W^, and H* is the 
charged Higgs boson. Note in the MSSM, due to the absence of S there only exist two CP- 
even Higgs bosons and one CP-odd Higgs boson, and consequently, U¥ and U^ are reduced 
to 2 x 2 matrices parameterized by the mixing angles o and 5 respectively. In our study, 
we choose the input parameters in the Higgs sector as (tan 3, ma, u) for the MSSM, (A, x, 
tan B, Hert, Ma, Ax) for the NMSSM with m3, = us (An + #), and (A, tan B, Herp, Ax, 
mg, ma) for the nMSSM with m2, = sas (HA, + r€pM?). 

The Yukawa couplings of the neutral Higgs bosons to the top and bottom quarks are 
given by [10 
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with U , U^ defined in Eq.(8). Obviously, once UZ / cos 8 « 1 as discussed later, the width 


13, 


After diagonalizing the neutralino mass matrix in the nMSSM, its mass takes the form |14 
2uA?* (v2 + v?) tan B 
which implies that Y? must be lighter than about 60 GeV for u > 100 GeV (required by 


LYukawa = 


Uigbysba, (9) 


of h; — bb is to be suppressed. 
Note the properties of the lightest neutralino X? in the nMSSM are quite peculiar 


(10) 


chargino mass bound) and A < 0.7 (required by perturbativity). If X? acts as the dark matter 
candidate, a light CP-odd Higgs boson a is then preferred to accelerate X? annihilation to get 
the acceptable dark matter relic density 13} Detailed study indicates that myo < 37 GeV 
and for most cases, m, < 60 GeV, which implies the SM-like Higgs boson ^ may decay into 


KKL or aa so that Br(h — yy) is suppressed|13}. 


4 


III. NUMERICAL RESULTS AND DISCUSSIONS 


A. Description of calculations 


To compare the signal rate with the SM prediction, we define a normalized rate as 


RSUSY 


osusy (pp > h > yy)/osu(pp > h > yy) 


K 


T(h > gg) Br(h > yy)|/T (hsm — gg) Br(hsu > y7)| 


[F( > gg)F(h > yy)]/[F(hsu > g9) (hsm > v)] X Tiøolhsm)/LFialh) (11) 


where we used the narrow width approximation and the fact that at leading order the cross 


section of the parton process gg — h is correlated with the decay width of h — gg by 
ie in T À 
6(gg — h) = comh ð (s — mj) = z—T(h 9g) 5(8 — mh). (12) 


In SUSY, the hyy coupling arises mainly from the loops mediated by W-boson, charged Higgs 
boson, charginos and the third generation fermions and sfermions, and the hgg coupling only 
from the loops mediated by third generation quarks and squarks. Consequently, the widths 
of h > yy, gg are given by 
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represent pure SUSY contributions with m; and m,, being sfermion mass and chargino mass 


respectively. Noting the asymptotic behavior of A" in the limit 7; « 1[15 


1 4 
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one can easily learn that the effects of the third generation squarks on the hyy and hgg cou- 
plings drop quickly as the squarks becomes heavy, and that the charged Higgs contribution 


to hyy coupling is usually far smaller than the W-boson contribution. 
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In SUSY, the third generation squarks can also affect the masses and the couplings of 
the CP-even Higgs bosons by radiative corrections, and such effects are maximized in the 
so-called “maximal mixing” (m7) scenario defined as X, = 2Msusy and A, = Ay, in the 
on-shell scheme b3. where X, — A, — u/tan with A, denoting the trilinear couplings 


of the top squarks and Msysy standing for the common soft breaking mass for the third 


generation squarks, i.e., Mg, = Mu, = Mp, = Msusy. Since the corrections are vital for 


our results, we will specially discuss them later. 
Different from previous studies in , we consider more constraints on the models, 


which are: 


(1) The constraints from the LEP-II direct search for neutral Higgs bosons in various 


possible channels. 


(2) The direct mass bounds on sparticles and Higgs boson from LEP and the Tevatron 


experiments |16}. 


(3) The LEP-I constraints on invisible Z decay: F(Z — X?2X9) < 1.76 MeV, and the 
LEP-II constraints on neutralino productions o(e*e^ — X9X?) < 107? pb (i > 1) and 


c(e*e- — X2X2) < 107 pb (i,j > 1) LT. 


(4) The indirect constraints from B-physics (such as b — sy) and from the precision 
electroweak observables such as pe, sin? m and My, or their combinations e;( = 
1,2,3) hs. We require €; to be compatible with the LEP/SLD data at 95% confidence 
level. We also require new physics prediction of Rẹ = I'(Z — bb) /I'(Z — hadrons) is 
within the 2c range of its experimental value. The latest results for Ry are Rọ” = 


0.21629 + 0.00066 and R7 = 0.21578 for m, = 173 GeV : 
(5) The constraints from Tevatron experiments on c(pp > h + X > 4y,2u27)|19]. 


(6) The constraints from the muon anomalous magnetic moment: af"? — gm ex (25.5 £ 


8.0) x 1077? [20]. We require the SUSY effects to explain a, at 2c level. 


(T) Dark matter constraints from the WMAP relic desity 0.0975 < QA? « 0.1213 bi. 
For each model we assume the lightest neutralino as the only component for the dark 


matter. 


As verified by numerous studies, these constraints show strong preference on the SUSY pa- 
rameters, e.g., the constraint (1) favors heavy top squarks with significant chiral mixing, 
while the constraint (6) favors large tan 8 for moderately heavy sleptons. Note that most 
of the constraints have been encoded in the program NMSSMTools b3l. which computes 
various Higgs decay rates up to one-loop level (the dominant one-loop and leading logarith- 
mic two-loop corrections to the Higgs masses and mixings are also included). We extend 
the code by adding more constraints in item (4) and further make it applicable to the 
nMSSM [13] (through some helpful discussions with the authors of the NMSSM Tools). 
Since the LHC is now testing the probability of the enhanced di-photon signal, we inves- 
tigate the situation where the signal rate can exceed its SM prediction. Eq. (QI) indicates 
two mechanisms in doing this. One is to enhance the hyy coupling or the hgg coupling. 
However, as indicated by our numerical results, this mechanism can only enhance the cou- 
plings by a factor up to 1.3 and 1.1, respectively. The reason is that the relevant SUSY 
parameters, such as tan 8 and the third generation squark masses, have been limited by the 
constraints. The other mechanism, which proves to be capable in enhancing RUSY by a 
factor up to 5, is to suppress the width of h — bb to enhance the branching ratio of h > ^y. 
To understand this, let's look at the expression of '(h; — bb), which, after including the 
important SUSY correction to bottom quark mass Ap, is given by|25 
Ui ) ( t U#/UE cot zx 
cos 8 1+A, 


where the first factor comes from the bottom Yukawa coupling in Eq.(9) and the second 


T'(h; — bb) x ( (17) 


factor arises from transforming the Higgs fields from weak basis to mass eigenstates in the low 
energy effective Lagrangian. Obviously, once U;2/ cos 8 « 1 and/or UJ /U3 cot BA > —1, 
T(h; — bb) will be greatly suppressed. In the following, we take the MSSM as an example 
to discuss how to satisfy the conditions. 

In the MSSM, Eq.(17]) may be rewritten as 


sina V^ f1— cota cot BA, V? 
cos 8 1+ A, 


where a is the mixing angle of the two CP-even Higgs boson obtained by diagonalizing the 


['(h — bb) œ ( (18) 


corresponding mass matrix M?,, and A, is given by 


Aj = ASCP 4 ABEW 


2o,mg Y? A: 
utan B (Stm, m, ma) + 562 


Tong mgs) +); 


T 


with the function J defined by 


1 nis 9. anu D 2921 c 
Habd = (Sa — eai — d a'b log qz + b'c log -zz + c'a UB . (19) 


m2, sin? B + m2, cos? 8 + Ay, — (m?, + m2) sin B cos B + A 
W a B Z B u (mi z) sin B cos B 12 (20) 


— (m2, + m2.) sin B cos B + A43. m2, cos? 8 + m2, sin? B + Ago 
A Z A Z 


where A;; (i, j = 1,2) denote the important radiative corrections with their leading contribu- 
Mi Mig 
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tions proportional to T In , one can numerically check that without A;;, sina//cosB 
Ww 


is always larger than unity for m B > 7 as required by muon anomalous momentum. So 
to suppress sino or equivalently the off-diagonal entry of the mass matrix the radiative 
correction must be present, and a positive large A1» along with a light CP-odd Higgs boson 
is efficiency in doing this. Meanwhile, given cot acot 6 ~ 1, A, must be around unity to 
satisfy cota cot GA, — 1, which requires large utan 8. In summary, in order to suppress 
T(h; — bb), light A as well as large u tan f is favored for given sparticle spectrum. We note 
what we are discussing is actually the so-called ‘small agg scenario’ of the MSSM [23]. 

From Eqs. QPO) one can infer that, in the heavy sparticle limit, the effective hyy and 
hgg couplings approach to their SM predictions and R^U9Y is determined by I'(h; — bb) 
or more generally by the total width l'i;(h); while in a general case, the contribution from 
the sparticle-loops to the couplings may interfere constructively or destructively with its 
corresponding SM contribution, and the size R?^"** then depends on the competition of 
D(h > gg) (h > yy) with 'i4(h). We checked that this conclusion is also applicable to the 
NMSSM and the nMSSM. 


B. Results for the MSSM in a general scenario 


To study R?"5Y quantatively we scan over the MSSM parameters under the constraints 
(1-7) and calculate the di-photon signal rate for the samples surviving the constraints. Since 
the first two generation squarks have little effects on the di-photon signal rate, in our scan 
we fix their soft parameters at 1 TeV. As for sleptons, since it only affects significantly the 
muon anomalous magnetic moment a,, which can in turn limit the important parameter 
tan B, we assume all soft parameters in slepton sector to take a common value m; and treat 


m; as a free parameter. For simplicity, we also assume the grand unification relation for the 
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gaugino masses, 3M1/5o4, = M2/a2 = M3/a3 with a; being the fine structure constants of 


the different gauge groups. Our scan regions are 


1 < tan 8 < 60, 90 GeV € m4 <1 TeV, 


200 GeV < Msusy ( Mas My; Mp;) <1 TeV, 


—3 TeV € Arp € 3 TeV, 100 GeV < u, Mo, m; € 1 TeV. (21) 


In Fig[I]we display the surviving samples, showing the di-photon signal ratio R“°°™ defined 
in Eq. (I) and the Higgs decay branching ratio versus the mass of the SM-like Higgs boson. 
This figure shows that in the MSSM there exist some points where R is enhanced by a factor 
up to 1.5. Such an enhancement is mainly due to the suppression of the total width of h, 
or equivalently the enhancement of Br(h > yy), which is shown in the right frame of Fig [i] 
Note that we required u < 1 TeV in our scan. If we relax u < 2 TeV in the scan, we find 
that RM?5M can be as large as 4. We checked that those samples giving R > 1 actually 
correspond to the ‘small og scenario’ discussed in , which is characterized by a large 


p tan f and |sin aesf/ cos B| < 1. 
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FIG. 1: The scatter plots of the surviving samples, showing the di-photon signal ratio RM55M 
defined in Eq.(11) and the Higgs decay branching ratio versus the mass of the SM-like Higgs 


boson. 


Fig [also shows that for most of the samples, the rate of the di-photon signal is suppressed 
relative to its SM prediction. These samples are usually characterized by an enhanced Abb 


coupling and a reduced hgg coupling (the change of the hyy coupling is usually negligible). 
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We checked that for R?"U5Y > 0.6 the effect of the reduced hgg coupling may be dominant, 
while for R5U9Y < 0.5 the effect of the enhanced hbb coupling is always dominant. We 
emphasize that for the samples with R5U5Y < 0.5, A must be relatively light (m4 < 300GeV) 
to ensure that the properties of h significantly deviate from the SM Higgs boson |26]. 

We note that current experiments can not rule out a light A with 110 GeV < m4 < 
140 GeV in the MSSM [27]. In this case, both A and H (the heavier CP-even Higgs boson) 
give rise to the di-photon signals similar to the SM-like Higgs boson h. However, the rates of 
these signals from A and H can not be large. This is because for 110 GeV < ma < 140 GeV, 
tan must be larger than 7 as required by the constraints (particularly by a,,) , which 
implies cosa > 0.8 from the tree-level relation tan 2a = tan 20 mate Since the Abb and 
Hbb couplings are proportional to tan 3 and cosa/cos 8 respectively, the branching ratios 


of A, H — yy are suppressed and so are their induced di-photon signals at the LHC [28]. 


C. Results for different models in the mj" scenario 


Since the NMSSM and the nMSSM have more free parameters than the MSSM, it is 
difficult to perform a general analysis of the signal rate. However, considering our aim is 
to show the differences of these three models, we examine the signal in the so-called m7*” 
scenario described in Sec. In this scenario, under the constraints (1-7) we scan over 


the following parameter ranges: 


90 GeV < my < 1 TeV,1 < tan 8 < 60, 100 GeV < u, Mo, m; € 1 TeV, 
100 GeV < Msusy ( Mo; My; Mp;) <1 TeV, (22) 


for the MSSM, 


0<A,K € 0.7, 90 GeV € ma € 1 TeV, 


100 GeV < Msusy ( Moa; Mu, Mp) < 1 TeV, 
1<tan 6 < 60, |A,| < 1 TeV, 100 GeV < u, M», m; < 1 TeV, (23) 


for the NMSSM, and 


0.01 < A € 0.7, 100 GeV < ma, u, Ma < 1000 GeV, 


100 GeV < Msusy ( Mq, Mu, Mp) < 1 TeV, 
1<tan8<60, — 1 TeV < A, < 1 TeV, 0< mg < 200 GeV, (24) 
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for the nMSSM with the soft parameters to be 100 GeV for the (v,,, i) sector in order to 
satisfy the a, constraint|13]. For other insensitive parameters we adopt the same assumption 


as in the last section. 
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FIG. 2: The scatter plots of the surviving samples in the m7'^* scenario of the MSSM, NMSSM 
and nMSSM, showing the di-photon signal ratio defined in Eq.(11) versus the mass of the SM-like 


Higgs boson. 
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FIG. 3: Same as Fig.2, but projected on different planes for the NMSSM. Here only the samples 


satisfying RNMSSM > 1 are plotted. 


Max 


In Fig [2] we show the di-photon signal rates in the m7'^* scenario for three models. This 
figure shows that in the nMSSM the signal is always suppressed. In the MSSM the signal 


is mostly suppressed, but in a tiny part of the parameter space the signal can be slightly 
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FIG. 4: Same as Fig.2, but showing the branching ratios of the SM-like Higgs decay to yy (‘e’, 


red), to aa (‘x’, black), and to X9X9 (‘o’, sky blue) in the nMSSM. 


enhanced. In the NMSSM, however, the signal can be enhanced in a sizable part of the 
parameter space (the enhancement factor can be as large as 7). In order to figure out the 
reason for such a large enhancement in the NMSSM, we concentrate on the samples with 
R > 1 and study the ratio o ^M55M /55M (pp — h) and Br(h — bb, yy). Our results are 
shown in Fig[3] which indicates that the production rate can be enhanced maximally by a 
factor of 1.25, while Br(h — yy) can be enhanced to 2 x 107? once Br(h — bb) is suppressed 
to several percent. This conclusion justifies our previous analysis about the mechanisms to 
enhance the signal. 

Quite surprisingly, we found that in the NMSSM the samples with R > 1 are unneces- 
sarily accompanied by a large u. The fundamental reason is that in SUSY the hbb coupling 
is determined by the Hg component of h, and in the NMSSM, due to the presence of sin- 
glet field component in h, the Abb coupling can be suppressed more efficiently than in the 
MSSM. We also noticed that once Br(h — bb) is suppressed, Br(h > VV*) (V = W, Z) 
may also get enhanced, which should be limited by the combined search for Higgs boson at 
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the Tevatron B. We checked that for Br(h > yy) ~ 107, Br(h + VV*) can be enhanced 
by a factor of 4 relative to its SM prediction. 

For the samples with a suppressed di-photon rate in the NMSSM and the nMSSM, we 
find that l'i4,(h) is usually enhanced (due to the enhanced Abb and/or the open-up of new 
decay modes ) and the hgg coupling is reduced. We checked that for R^U5Y < 0.5 the 
former effect is dominant. We note that in the nMSSM R?U*Y is usually small, which is 
mainly due to the open-up of new decay modes of h, such as h + X9X9 (i = 1,2) or h > aa 
with their rates shown in FigH] We emphasize that this feature comes from the peculiarity 
of x? in the nMSSM (see Eq.(10)) and should keep valid regardless our choice of the soft 
parameters in the squark sector. We numerically checked this point by a more general scan 
than Eq.(24). We also note that for nearly all the samples in the NMSSM with m, > 120GeV 
we have RU5Y < 1, and for all the samples in the nMSSM with m; > 125GeV we have 
RSUSY < 0.14. We owe this to the constraints we considered, which severely constrained 
the enhancement of the branching ratio of h > yy (see Fig] and Fig). 

We also studied the di-photon signal rate in the ‘no-mixing’ scenario defined as A, = Ap 
and X, — 0. However, we found it is difficult for this scenario to satisfy the constraints if 
Msusy « 1 TeV, especially we did not find any surviving samples for the MSSM. Since the 
di-photon signal for the surviving samples in the NMSSM and the nMSSM do not exhibit 
new characteristics, we do not present the results here. 

So, we see that in low energy SUSY, depending on the models, the di-photon signal rate 
at the LHC may be significantly suppressed or enhanced relative to the SM prediction. With 
2fb ! integrated luminosity at the running LHC, the di-photon signal can allow for a test 
of the SM and a probe of the low energy SUSY models. For example, if the di-photon signal 
rate is found to be not smaller than the SM prediction, then the nMSSM will be immediately 
excluded (note that in this case the universal extra dimension and the little Higgs theory 


will also be ruled out because they suppressed the diphoton signal rate |6, |29]). 


IV. CONCLUSION 


We focused on the di-photon Higgs signal gg > h — yy for the SM-like Higgs boson at 
the LHC and performed a comparative study for three SUSY models: the MSSM, NMSSM 


and nMSSM. Considering various collider and cosmological constraints, we scanned over the 
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parameter space and obtained the following observation in the allowed parameter space: (i) 
In the nMSSM the signal rate is always suppressed; (ii) In the MSSM the signal rate is 
suppressed in most cases, but in a tiny corner of the parameter space it can be enhanced 
(maximally by a factor of 2); (iii) In the NMSSM the signal rate can be suppressed or 
enhanced depending on the parameter space, and the enhancement factor can be as large 
as T. 

Note added: After we finished the manuscript, we noticed a preliminary result from the 
ATLAS collaboration , which excluded R > 4.2 (R is defined in Eq[11] for m; ~ 115GeV. 
This means that in the middle panel of Fig.2 the samples above R ~ 4.2 for the NMSSM 


will be excluded. 
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